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Abstract 

We discuss a realisation of the pentaquark structure proposed by Jaffe and Wilczek within a 
simple quark model with colour-spin contact interactions and coloured harmonic confinement, 
which accurately describes the A — N splitting. In this model spatially compact diquarks are 
formed in the pentaquark but no such compact object exists in the nucleon. The colour- 
spin attraction brings the Jaffe- Wilczek-like state down to a low mass, compatible with the 
experimental observation and below that of the naive ground state with all S-waves. We 
find, however, that although these trends are maintained, the extreme effects observed do not 
survive the required "smearing" of the delta function contact interaction. We also demonstrate 
the weakness of the "schematic" approximation when applied to a system containing a P- 
wave. An estimate of the anti-charmed pentaquark mass is made which is in line with the 
Jaffe-Wilczek prediction and significantly less than the value reported by the HI collaboration. 
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1 Introduction 



In the wake of the possible discovery of baryon with positive strangeness (the O) [1], phenomeno- 
logical models have arisen designed to explain its low mass (~ 1540MeV) and narrow width 
(< lOMeV). For an early review see [2]. 

One such model is due to Jaffe and Wilczek [3] . They propose diquark correlations of low mass 
in a relative P-wave, giving the state overall positive parity in line with the Chiral Soliton Model 
prediction [4]. 

The Jaffe- Wilczek model has two diquarks each of spin-0, flavour 3f and colour 3 C , which by 
Bose statistics must couple symmetrically. To give an overall 10^ the diquarks must couple to 
a 6f , i.e. symmetrically. An overall colour singlet requires antisymmetric diquark coupling to a 
3 C and hence to be overall symmetric we are forced to introduce a P-wave between the diquarks. 
A P-wave appears in many quark model treatments of pentaquark structure, but, unlike in the 
Jaffe- Wilczek model, it is usually rather ad hoc. 

This model has attracted much attention, and is particularly appealing as the same diquark 
correlations can be used to explain the enigmatic scalar mesons below lGeV as being diquark- 
antidiquark states [5,6]. 

Our point of departure is to discuss the dynamics which gives rise to diquarks. This an open 
question which must eventually be answered directly from QCD, but we will show that such 
correlations can appear even within simple quark potential models of the type often used to 
describe the conventional hadron spectrum. 

2 Schematic Approximation, Quark Potential Model and a 
Jaffe-Wilczek-like state 

Previous investigations into pentaquark structure have often made use of the "Schematic" ap- 
proximation to colour-spin forces [7]. This approximation discards spatial dependence, having an 

interaction potential Vg CH — — C^2 i . ^ • -^-§i.Sj, where C is a constant. For example in the 
total spin-0 channel this would give, 



where e.g. 1 ® 1 — > indicates that 1 & 2 and 3 & 4 are each coupled to spin-1 and then the two 
pairs coupled to total spin-0. 

We will go beyond this approximation by allowing there to be a non-trivial spatial dependence. 
In particular we shall see that the schematic approximation, while capturing all the physics for 
spatially symmetric states, is not sufficiently versatile to accurately describe states containing a 
P-wave. This will be demonstrated by computing the equivalent to Q in a model with non-trivial 
spatial dependence. 

We now introduce the model which, in essence, is a standard quark potential model suitable for 
describing the light baryon and meson spectrum. 

As a binding potential between quarks we take the coloured harmonic oscillator, 
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Figure 1: Five-body co-ordinate set in the q 4 



q rest frame. 



This is chosen mainly for ease of calculation but has been used in the past as an approximation 
to the more phenomenologically justified coloured Coulomb + linear potential. 

In addition to this we introduce a colour-spin contact interaction, 



It is this potential that will act to bind two quarks tightly into an effective light diquark. 

The fundamental mechanism binding the diquarks in the Jaffe-Wilczek picture is unlikely to 
be anything so simple as a colour-spin contact interaction, and indeed it may not be possible to 
consider diquarks in terms of constituent quarks with non-relativistic interactions at all. Diquarks 
may appear as degrees-of-freedom in the same way as constituent quarks do during chiral symmetry 
breaking. Despite this, colour-spin contact interactions have a considerable phenomenological 
pedigree and may well be able to simulate quite effectively the true QCD interaction responsible 
for the "fine" structure in hadrons and as such we will investigate their effect on the pentaquark 
system. The motivation we have for a delta-function interaction is that of [8], the Breit-Fermi 
reduction of one-gluon exchange between quarks, but delta-function interactions have also been 
used to model instanton effects in light hadrons (see for example [9] and references therein). 

Working in the q A q rest frame, we define the internal variables, r, A, p, v 7 by 



V"(q iqj ) = ■ $f Si.Sj 5(rij). 
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which are clearly well suited to describing a Jaffc-Wilczek-like configuration if the diquarks are 



spatially correlated (and see figure With these variables the kinetic energy is 2 
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and we can write the Hamiltonian H(q A q) = T(q A q) + V(q A ) + V{q) + V a {q A ) + V a (q). 

We can approximately solve for an eigenstate of this Hamiltonian using the variational method. 
First introduce a trial wavcfunction containing a number of variational parameters and minimize 
the expectation value of the Hamiltonian with respect to them. The trial wavcfunction has spatial 
form 

V m = ^[/3Ar lm (A)e-^ A 'V 2 ][r 00 (p) e ^^V2 nroo( p )e -7^V2 H y oo(f)e -aV 2 /2 ]j (6) 

which is a P-wave in the A variable and an S-wave in the others in the harmonic oscillator 
approximation. The flavour-spin-colour structure is chosen to be |3o,3 )g ® |3 c ,3 c ) 3c , where 
this is a shorthand for quarks 1 & 2 coupled to a 3 of flavour, spin and colour 3; quarks 3 & 
4 are coupled identically to 1 & 2 and then the pairs {12}, {34} coupled to a 6 of flavour, spin 
and a 3 of colour. This is the flavour-spin-colour structure of the Jaffe-Wilczek correlation. The 
antiquark couples trivially to give an overall colour singlet. 

As stated this wavefunction (flavour-spin-colour-spatial) is only antisymmetric under exchange 
of labels 1 «-> 2 or 3 «-> 4, but not for example 1 <-> 3. This would seem to violate the generalized 
Pauli principle. In fact we need not go to the trouble of antisymmetrising two particles if their 
wavepackets have very limited overlap - we don't have to worry about electrons on the moon when 
we solve the Schrodinger equation for a hydrogen atom on Earth nor, more pertinently, do we 
usually antisymmetrise the quarks in different nucleons when we study the deuteron. Here we are 
proposing that the dynamics in the Hamiltonian given earlier will be such that the wavefunctions 
of quarks in different diquarks have very limited overlap such that we don't have to antisymmetrise 
them. We will show that this assumption is consistent when we variationally solve the Hamiltonian 
with this trial wavefunction. 

With the |3 c ,3 c )3 c colour state the expectation values of the potentials are 
(3 c ,3 c |y(g 4 )|3 c ,3 c ) = b - a {p 2 + v 2 ) + ^a\ 2 
(3 c ,3 c |F(g)|3 c ,3 c ) = ia(p 2 + v 2 ) + ^aA 2 + ^ar 2 
(3 c ,3 c \V°(q A )\3 c ,3 c ) = |/i(5i.5 2 <5(ri 2 ) + S 3 .S 4 5(r 34 )) 

+ ~h(Si.S 3 8(r 13) + Si.S^rU) + S 2 -S 3 5(f 2 3) + ^.S^Sfaij) 
(3 c ,3 c \V°(q)\3 c ,3 c ) = h % 4 ).S ? (g)<5(r l9 -). (7) 

2 11 is the q 4 , q reduced mass, m g 4m 
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With the spatial wavefunction © we have 
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The spin structure of (3 C , 3 c |V rCT (g 4 )|3 c , 3 C 
spin-0 correlation, 
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guarantees an attractive potential only for the spin-0, 
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Compare this with the result of the schematic approximation, eqn([Q). We see that by suitable 
choice of C, the schematic approximation can duplicate the <g> result, but that unless there 
is accidental cancellation it would be unable to describe the 1 ® 1 — > result. The origin of the 
non-zero value for 1 <g> 1 — > is the different expectation values of the delta function with a P-wave 
and without (see equation©). We propose that this is a general problem with the schematic 
approximation when applied to states with non-trivial spatial dependence and that by its use one 
can miss significant physics. 

Returning to our study of the quark potential model we see that (3 C , 3 c |F°'(g)|3c, 3 C ) is particu- 
larly simple. The spatial and colour dependence is identical for all quarks and factors out leaving 
a sum of the quark spins, which in the £g> — > channel is zero, hence the anti-quark does not 
change the hyperfine energy of the system in the Jaffe-Wilczek correlation. 

The expectation value of the Hamiltonian is thus, 

\3a 2 2a] T5/3 2 5a 1 „ T37 2 3al / 7 \ 3 
{ JW\H\JW)=[— + ^\ + [£ + w \ + 2[-^ + - 2 \-h^ . (10) 

The minimum is found to be at a = (S/ia/3) 1 / 4 , (3 = (2ma) 1 / 4 and 7 satisfying ^ 2 ^ /2 7 5 — ^;7 4 + 
12a = 0. We can set the parameters m, a, h using conventional baryon spectroscopy and along the 
way demonstrate that this model can perfectly well describe the A — N splitting. 

a can be set approximately using the S- wave P-wave splitting (ujp) of roughly 500MeV for the 

non-strange baryons, since in the harmonic oscillator ujp — y/ia/m. This gives a ~ 4^330 m2v ' 
where the reason for this unusual presentation will become clear later. If we consider one-gluon- 
exchange to be the origin of the contact term [8, 10], then h = ^-^z, where as ~ 0.75 is usual. 
The light quark mass takes its conventional value m = 330MeV. 

The expectation value of the Hamiltonian for the nucleon between 5*-wave Gaussian trial wave- 
functions (with oscillator parameter a p ) is then 



which is minimised by a p = 440MeV, with colour-spin hyperfine energy — 150MeV. This corre- 
sponds to a A — TV splitting of 300MeV, in good agreement with data [11] 3 . My approach differs 
from that usually taken which considers the hyperfine term only in first order of perturbation the- 
ory. By using a variational ansatz we are allowing the hyperfine term to modify both the energy 
and the wavefunction (a p would have been 405MeV had the hyperfine term been neglected - hence 
the rather unusual form of a presented earlier which makes this obvious) . 

Returning to the Jaffe-Wilczek state, using these parameter values we find a — 370MeV, (3 = 
340MeV and 7 = 520MeV. Had we neglected the hyperfine term, 7 would have been 405MeV. The 
hyperfine term is reducing the mean distance between quarks 1 & 2 and quarks 3 & 4 relative to 
the {12}{34} distance, which is exactly what one demands of a spatial diquark-diquark state and 
what we need for the unsymmetrised ansatz to be justified. Specifically with these numbers the 
diquarks have mean radius ~ 1/3 fm and are separated by an average distance ~ 1 fm. 

With 7 = 520MeV, the hyperfine energy is — 510MeV, much larger than the cost of exciting the 
P-wave ~ 350MeV), so that the Jaffe-Wilczek state is considerably lighter than one would 
have naively expected. Simply adding together quark masses and including a P-wave energy of 
~ 350 MeV would suggest a mass over 2 GeV. Consider the difference m(0) — m(N) which in this 
model (where the Hamiltonian does not include the quark rest masses) will be 

( JW\H\ JW) +m + m s - (N\H\N). 

This is found to be (with m s = 450MeV), ~ 590MeV and hence m(8) <~ 1530MeV. This agreement 
is fortuitous - we have, for example, neglected tensor interactions which are non-zero due to the 
P-wave character and which will change this prediction. 

We have observed something rather interesting - A Jaffe-Wilczek-like 9 state of low mass with 
compact diquarks has emerged from a model which also gives an excellent description of the A — N 
splitting. From conventional quark model analysis one would not have expected this - the usual 
colour-spin arguments as applied to the nucleon suggest a spin-zero "diquark" of mass <~ 600 MeV 
which if used naively in the Jaffe-Wilczek scheme would hugely overpredict the 6 mass. What we 
have demonstrated here is that the diquarks in the O are not like the "diquark" in the nucleon. 
The quarks in the nucleon all overlap spatially and correct antisymmetrisation must be carried 
out. This does not allow for the kind of spatially distinct diquark found for the 9, instead the 
correlation can only be in spin and flavour. 

3 Flavour-Spin 210 "ground state" 

The Jaffe-Wilczek correlation would not be our first guess for the ground state of the q 4 q system. 
Without colour-spin interactions a state with S-waves between all quarks would be expected to 
have lower energy. Jaffc & Wilczek suggest that the colour-spin interactions would force such a 
state to a higher energy than their correlation and it is this to which we now turn in this simple 
model. 

The totally antisymmetric q 4 state with symmetric spatial part is in a 210 of flavour-spin and 

3 of course, the parameter ag has been chosen to give this good agreement 
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is coupled to total spin S(q 4 ) = 1 [12]. It has explicit form 

-J= ||6i,6i) ® |3c,3 c ) + ^|3i,3 ) - ||6 ,6i>^ ® |6 C ,3 C ) + ^|3 ,3i) - ^|6i,6 )^ ® |3 C ,6 C 

(12) 

where we use the "diquark" notation without implying that diquarks are dynamically generated. 
The spatial form (for q 4 ) is 

V = Ar[yoo(A)e-^^/2 ][ y oo(/5)e -^, 2 /2 ][ y oo( ^ )e -^^/2 ]5 (13) 

whose overall symmetry is exposed by expressing the exponent in terms of the n as — 2(3 2 J2i>j — 
rS) 2 - 

In this flavour-spin-colour state the potential V(q 4 ) is -^(A 2 + p 2 + v 2 ). A straightforward but 
somewhat lengthy computation yields the expectation of V a (q 4 ) in this state. This is simplified 

slightly by the expectation of the delta function, (S(fij)) = (-J^j i being i, j independent. We 
find 

which is repulsive, as anticipated. Comparing with the q 4 part of the Hamiltonian in the Jaffe- 
Wilczek case we find that the 210^5 is around 60MeV heavier. We have not considered the 
colour-spin interaction between q and the anti-quark which may be non-zero due to the S(q 4 ) = 1 
character of the 210_fs state. This could raise the mass of the 210_fs even further. 

Thus in this simple model with colour-spin contact interactions a Jaffe-Wilczek-like state can 
be consistently defined which is lighter than the naive ground state. The diquark size is small on 
the scale of their separation, which was assumed in the spin-orbit analysis of [13]. 
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4 Smearing the delta function 

One should worry about the validity of using a delta function interaction in a Hamiltonian. Such 
an object is too singular at the origin to use the normal boundary condition u(0) = to quantise 
the energy levels; we have ignored this problem by using a variational ansatz which satisfies the 
usual boundary condition. We still have a problem if we consider the origin of this term - it came 
at order (v/c) 2 in a non-relativistic reduction of the one-gluon-exchange process, hence in using 
this term we are implicitly assuming that the momenta of the quarks are much lower than their 
mass. At very small interquark separations the corresponding momentum scale is rather large 
and the non-relativistic approximation breaks down. As such we should "smear out" the delta 
function on the scale of the quark mass. A suitable modification is 

^-^jexph^/rg] (15) 

with ro ~ m~ . This type of smearing has been used previously in the literature when modeling 
conventional hadron spectra [9, 14]. The following smearing scales were found to be phenomeno- 
logically satisfactory: 
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[9] r ~ l/(1280MeV) ~ l/(4m„); 
[14] r ~ l/(1870MeV) ~ l/(6m u ). 
We will consider here the effect of smearing on the results we reported in the previous section. 
The change is (0\5{r)\0) = {0/^f {0/^f{rl0 2 + 1)" 3 / 2 . With r = l/(ran u ), the following 
is obtained for the nucleon, Jaffe-Wilczek state and the 210ps 
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The delta-function is recovered in the limit n — > oo. So while the nucleon and 210fs state are 
not strongly affected by the smearing for the phenomenological values n ~ 4 — + 6, the large effects 
felt by the Jaffe-Wilczek state are diluted considerably. Since the P-wave excitation energy is 
independent of n, the dilution is such that the Jaffe-Wilczek state is heavier than the 210fs- 

That the conclusions arrived at in the previous sections are cut-off dependent is disappoint- 
ing, but at least the trends remain. The Jaffe-Wilczek-like state still undergoes a considerable 
downward shift and spatially localised diquarks still seem to be formed. 

5 Extensions 

As mentioned in the introduction, an attractive feature of diquark correlations is their supposed 
ability to describe the light scalar mesons. A simple extension to this work would consider qqqq 
states; however, without a P-wave to separate the diquark from the anti-diquark it seems unlikely 
that spatially distinct diquarks of the type found earlier will emerge. Of course such a study would 
share many similarities with the classic analysis (in the bag model) of Jaffe [5] . 

Another interesting state to consider is the ududud state with two P-waves. This has the same 
quantum numbers as the deuteron. If three diquarks are formed with P-waves between them we 
would hope that it has mass greater than the measured mass of the deuteron which we know to 
be well described as a bound-state of a proton and a neutron with small overlap. Consideration 
of such a state has been advocated in [13] as a test on any model proposed to give a light 6. 

This model, extended to allow unequal quark masses, could predict the masses of the other 
members of the pentaquark 10® 8, being an explicit realisation of SU(3)f breaking. In particular 
this would test the phenomenological Jaffe-Wilczek Hamiltonian H s = Mq + (n s + n s )m s + n s a 
and its predictions for the other pentaquark states. 

If confirmed, the anti-charmed pentaquark observed by HI [15] with a mass of 3.1 GeV would 
have its mass significantly undcrpredicted by this model, as it is in the original Jaffc-Wilczck 
paper [3]. Setting the charm quark mass using the experimental A c — A mass difference and 
computing the new Hamiltonian expectation, we find a mass of around 2.8 GeV for the equivalent 
anti-charmed Jaffe-Wilczek-like state, which is within 100 MeV of the Jaffe-Wilczek prediction 
and very close to the DN threshold. If the magnitude of spin-orbit splitting calculated in [13] is 
correct, the possibility that the observed state is the 3/2 + state is unlikely. One possibility would 
be that the lightest anti-charmed pentaquark pair (l/2 + ,3/2+) is still to be found and that the 
HI observation is an excited state, such as the vector diquark excitation. The state 1 <g> 1 — > 
in the model presented in this paper is about 400MeV heavier than the Jaffe-Wilczek-like state. 
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This would be a little heavy for the HI candidate but probably within model errors; alternatively 
the state with one vector diquark and one scalar diquark will be somewhat lighter and might be 
a possibility, but only if the HI state is the I z = part of an isovector. 
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